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Abstract. A simple algebraic approach to the study of multipartite entanglement 
for pure states is introduced together with a class of suitable functionals able to 
detect entanglement. On this basis, some known results are reproduced. Indeed, 
by investigating the properties of the introduced functionals, it is shown that a subset 
of such class is strictly connected to the purity. Moreover, a direct and basic solution 
to the problem of the simultaneous maximization of three appropriate functionals for 
three-qubit states is provided, confirming that the simultaneous maximization of the 
entanglement for all possible bipartitions is compatible only with the structure of 
GHZ-states. 
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1. Introduction 

Entanglement plays an important role in applications of quantum mechanics to 
nanotechnologies, especially in the field of quantum information and communication. 
The concept of entanglement is mathematical, since it corresponds, in the case of pure 
states, to the impossibility of writing the quantum state describing a compound physical 
system as a simple product of states related to the single subsystems. Nevertheless, this 
mathematical property of the state has interesting physical implications, such as the 
presence of non classical correlations between quantum systems, even when they are 
quite far from each other (non locality). 

In spite of the strong interest and deep studies developed on this subject, detection 
and classification of entanglement in multipartite systems are unsolved problems up 
to date. Indeed, it has not yet been defined a universal quantity able to measure the 
entanglement level to be associated to a generic pure or even mixed state. (The concept 
of non entangled state is translated into factorizability for pure states and separability 
- simple or generalized — for mixed states, though the latter concept includes the 
former as a special case.) Besides, there are properties that any measurement of 
entanglement should possess. In particular, it must be discriminant (vanishing iff 
the state is separable), not increasing under LOCC (local operations and classical 
communication) and convex (it must not increase when two or more states are combined 
in a mixture) pp. 

Though a general answer to measurement and detection of entanglement is absent, 
if the system is composed by a few of low dimensional subsystems it is possible to 
provide conditions for the presence of entanglement, but rarely these conditions are 
both necessary and sufficient. In 1996, Peres [2] introduced a sufficient condition for 
the detection of bipartite entanglement that can be written in terms of a functional 
addressed as negativity, subsequently studied by Horodecki et al, who have proven that 
this condition is necessary and sufficient for low-dimensional systems [31 0]. Another 
important functional is the concurrence introduced by Wootters in [5j [6] , which is a true 
measure of entanglement, both for pure and mixed states, but unfortunately applicable 
only to systems which are couples of two level systems (i.e., a couple of qubits). 

The study of multipartite entanglement is more complex [7J EJ |9], not only from 
the computational point of view, but even at a conceptual level, since for instance it is 
neither immediate nor intuitive to understand what is a state with a maximum level of 
multipartite entanglement. In 2008, Sabin et al, [10], tried to extend the negativity 
to the tripartite case, succeeding in finding a functional that detects the tripartite 
entanglement when applied to pure states. However, this functional gives only a clue 
about the presence of entanglement in mixed states, indeed in this case the condition 
given is neither necessary nor sufficient. Nevertheless, its effectiveness has been shown 
in the study of simple physical systems [11]. Another interesting quantity, introduced 
for detecting entanglement of pure states of three qubits, is the three-tangle [12], which 
is based on the concurrence and whose validity has been criticized [13]. Recently, a 
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classification of entangled states for three-qubit states has been given on the basis of 
some invariant quantities [H]. Moreover, attempts to apply well known results of algebra 
or geometry to study multipartite entanglement have been made. In [15], for instance, 
Makela et al associated to every pure state of N-qubits a polynomial, characterizing 
completely the state, capable not only to detect factorizability of the state but even the 
number of separable qubits. Instead, in [16], Miyake used the hyperdeterminants and 
the theory of Segre variety to give a classification of the entanglement of pure states in 
tripartite systems. A very useful tool to study multipartite entanglement in mixed states 
is given by the necessary condition for separability expressed by Huber et al [17] , which 
exploits a correlation function defined through replicas of the states under scrutiny and 
operators of partial swapping, and which has been exploited to reveal thermal tripartite 
correlations in spin-star systems |18j . 

Another important quantity that has been used to reveal entanglement in several 
cases (bipartite and multipartite systems), is the purity (strictly connected to the linear 
entropy [IH]) of the reduced density matrix, which in passing is strictly connected with 
the concept of mixedness [20] : the more two systems are entangled, the less pure is the 
reduced state describing each one of the two systems. 

In this paper we reproduce some known results about entanglement detection 
by introducing a simple algebraic approach to establish whether a pure state of a 
multipartite system is entangled with respect to a given bipartition. Through this 
analysis it is possible to determine whether a given pure state is completely separable, 
separable or totally entangled. In the last case one can infer the presence of genuine 
multipartite entanglement. On the basis of this approach, we naturally introduce a class 
of functionals which includes, among others, quantities traceable back to the purity and 
the linear entropy, strictly related to the so called concurrence vectors, presented in [21] 
by Akhtarshenas and studied by Mintert et al in [22] [23] . 

The paper is organized as follows. In the next section we describe our approach 
to the study of multipartite entanglement and give the definition of a relevant class 
functionals, proving the ability of such quantities to detect entanglement in pure states. 
In section [3] we show that such quantities are strictly connected with purity and 
mixedness. In section H] we show, by exploiting our simple and algebraic approach, that 
in the case of three qubits, simultaneous maximization of the three relevant functionals 
is possible only for states equivalent (i.e., equal up to a local and unitary transformation) 
to the GHZ state. Finally, in the last section, we give some comments and concluding 
remarks. 

2. Factorizability conditions 

Each pure state of a bipartite system can be written in the form: 




(1) 



■j 
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where we take as a standard basis in the Hilbert space % = 7-L\ ®%2- This state 

is said factorizable if it can be written in the form: 

Ii>). (2) 

i 3 

This means that the state is factorizable iff i/3j, Vz,j. So, introducing the 

matrix of the probability amplitudes A = the state \<f>) is factorizable iff columns 
(respectively rows) span a 1-dimensional subspace, i.e. rank A = 1. 

In general, this means det A = and in the case of a two qubits system, in [21] it is 
proven that this condition is equivalent to the Peres-Horodecki criterion. The condition 
here expressed can be easily extended to investigate the presence of entanglement in 
every pure states of multipartite systems. We can summarize this result for bipartite 
systems with the following theorem. 

Theorem 1. A pure state of a bipartite m x n system is factorizable iff the rank of the 
corresponding matrix of the probability amplitudes is unity. 



Now we can extend this result to tripartite systems and we can generalize it to 
every multipartite system. Consider a m x n x p Hilbert space. Each pure state \<ft) 
living in it can be written in the form: 

m— 1 re— 1 p— 1 

\4>) = £££aijfclW> ( 3 ) 

i=0 j=0 k=0 

where {\ijk)} is the standard product basis for the Hilbert space related to the whole 
system. Again, the state is factorizable with respect, for instance, to the first component, 
iff: 

a ijk = aif3 jk Vi,j,k. (4) 

This time the matrix of the probability amplitudes is a cubic one and the 
separability condition of the first component can be traced back to the proportionality 
of the layers in the direction of %. For a better visualization of this statement we can 
write two layers in the z-direction: 



a 



a 



•Pip \ 



( a s /3u 

\ Ctsfinl 



a 



a 



■>Pnp J 



\ «r Ail • 

which are the r th and s th layers. 

Starting from these considerations, we can prove the following theorem. 

in a m x n x p Hilbert space is factorizable iff one of the 



Theorem 2. A pure state 
following is true: 

• CLijkO'i'j'k' — Clij'k'O'i'jk = 

• Q>ijk a i'j'k' — CLi'jk'O-ij'k = 



0, Vi,j,k,i',j' 7 kf; 
0, Vi,j,k,i',j',k'; 
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• Clijk a i'j'k' — a i'j'k a ijk> ~ , Vz, j, k, i' , j , k . 

Proof. For the sake of simplicity and without loss of generality, let us consider the 
factorizability condition for the first qubit. 

If the state \<f>) is factorizable with respect to such first qubit, using the 
factorizability condition in equation (J4]) we have: 

ctiPjkCti'Pj'k' - ctifij'k'Cti'Pjk = 0. (5) 

Vice versa, if we suppose aij^vyk 1 — dij'k'di'jk = 0, Vz, j, k,i',f, k' then, fixing i,j, k such 
that djjfc 7^ 0, it turns out that: 

O'i'j'k' — —Oij'k' — Cti'Pj'k' > Vi', (6) 

@>ijk 

where we take ay = and B-iy = au/fc/, which is the factorizability condition of the 

a ijk J J 

first component. □ 

Using this theorem we can define a class of functionals able to detect entanglement 
in each pure state of a tripartite system. 

Let / : C — > R be a positive function such that Vx G Cf(x) > and f(x) = iff 
x = 0, then we have the following: 

Theorem 3. A pure state \<p) in a m x n x p Hilbert space is factorizable iff at least 
one of the following is true: 



• m[ s) 




,k,i',j',k' f( a ijk^i'j'k' 


— aijik'Cii'jk) 


= 


. M 2 (/) 




,k,i',j',k' f( a ijkdi'j'k' 


— O'i'jk'Q'ij'k) 


= 


• mP 




,k,i',j',k' f( a ijkdi'j'k' 


— di'j'kdijk') 


= 



The proof of this result follows immediately from the theorem [2] and from the 
properties of /. 

( f) 

We can make some observation about this. Firstly, if Mj; = then the 
corresponding state is factorizable with respect to the k th component, so if two of these 
quantities are zero also the third has to be equal to zero (considering that, for pure 
states, if two of the three subsystems are separable, then the third one is separable as 
well). This statement can be easily proven considering the factorizability condition in 
equation (Tj0) for two of the three components. Secondly, if all the three are equal 
to zero than the state is completely factorizable — i.e., it is the product of three states 
of the three subsystems — and if each one of them is different from zero the state is 
genuinely tripartite entangled. 

Theorem [3] can be extended to the multipartite case considering all the possible 
bipartitions of the system. For instance, in the case of a quadripartite system we have 
seven conditions given by: 

• M i '■= Hi,j,k,i,i',j',k',i' f( a ijkiai'j'k'i' - aij'k'va-i'jki) = 0; 



An algebraic approach to the study of multipartite entanglement 



• 


M 2 (/) 




,l,i',j',k',l' f(, a ijklO>i'j 


k'V 


— aiijkn'O'ij'ki) 


= 


• 






lA'j'k'V f( a ijkl a i'j 


k'V 


— aiijiki'dijk'i) 


= 


• 


Mf 






k'V 


— 0'i'j l k'l a ijkl l ) 


= 


• 


Mg> 


^-si,j,k 


,l,i',j',k',l' f( a ijkl a i'j 


k'V 


— a>ijk'i ,a >i'j'ki) 


= 


• 


Mff 


^-si,j,k 


,l,i',j',k',l' f( a ijkl a i'j 


k'V 


— O'ij'kl'O'i'jk'l) 


= 


• 




^-si,j,k 


,l,i',j',k',l' f( a ijkl a i'j 


k'V 


— aiijki'O-ij'k'i) 


= 



where Mf: functionals are related to separation of the fc-th subsystem when the other 

( f) 

three are considered as a whole, while M^.j J, s are related to the bipartition obtained by 
considering the couple k — j as a whole and the other two parts as the second subsystem 
of the bipartition. 

3. Connection with the Purity 

In this section we prove that if we take as function / the squared modulus, the relevant 
functionals, simply denoted as are strictly connected with the purity of the relevant 
reduced density operators. To this end, let us first of all rewrite by considering 
only two indexes (this means that the two indexes, say i and j, will span multipartite 
subsystems) : 

M k = \aijai'j' — aij'ai'j\ 2 . (7) 

Then, by expanding the modulus, one gets: 

Mk = ^ y (aija*jai>j>a*/j l + aijia*j/ai>ja*ij — aijai'jia*j,a*>j — a*ja*ij f aijiai'j) 

= 2(l-p) (8) 



where 



V 



/j l >: l a ''i ,a '; i ' a ',\j - ^2 a *j a *'j' a ij' a i'j > (9) 



is the purity of the density operator associated to any of the two subsystems constituting 
the bipartition. This makes the functional Mf, essentially equal (up to the factor 2) to 
the linear entropy (mixedness) |20j . 

Therefore, the quantities Mk are invariant under local unitary transformations and 
are bounded as follows: < Mk < 2(D — 1)/D, where D is the dimension of the smaller 
among the two Hilbert subspaces associated to the bipartition considered. Moreover, 
maximization is reached when the whole-system state can be written as a superposition 
of D states with the same weights: 

l^) = ^ 1/2 Ee^|j)J$ J ),, (10) 

3 = 1 



An algebraic approach to the study of multipartite entanglement 



7 




Figure 1. Matrix of the probability amplitudes related to the W-state. The vertexes 
of the first layer correspond (from top-part left-side, clockwise) to the coefficients: aoooj 
aooi, aon, aoio; the vertexes of the second layer instead correspond to the coefficients: 

a 100> llOl; Glllj a no- 



where \j) k and refer to the two subsystems constituting the bipartition, and 

(j\l) = (Q j \*i) = 5 jl . 

3.1. Some Examples 

For a better understanding of the results given in the previous section we apply the 
previously proven criteria to some pure states. First of all, let us consider the W-state 
which, as well known, is genuinely tripartite entangled: 

|W0 = -^(|100> + |010> + |001». (11) 

The relevant cubic matrix of probability amplitudes is shown in figJU 
A rapid calculation gives: 

Mi = M 2 = M 3 = ~ (12) 
y 

and therefore, the state is genuinely tripartite entangled in accordance with theorem [3j 

As a second example, let us consider the GHZ-state: 

|G#Z> = i=(|000> + |lll>), (13) 
whose matrix of the probability amplitudes is shown in fig. [2j 
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Figure 2. Matrix of the probability amplitudes related to the GHZ-state. The vertex- 
coefficient correspondence is the same as in the previous figure. 



Once again a straightforward calculation leads to: 

Mi = M 2 = M 3 = 1. (14) 

By the way, we mention that in the following section we prove unity to be the 
maximum value that these quantities can reach for a three-qubit system. Moreover, it 
is worth noting that maximization of two does not imply maximization of the third, 
as shown by the following example. Indeed, consider the state: 



IVU = |ooo) + A |oii> , 

which is factorizable with respect to the first qubit, being Mi = 
factorizable with respect to the other two qubits, being Mi = M3 = 1. 



(15) 

0, but it is not 



Finally, let us consider a state for which it is not evident whether it is entangled or 



not: 



\<P) 



2V26 



(5 |000) + 3 |010) + 2 |001) + 4 |011) + 7 |101) + 1 1 11>), (16) 



Since it turns out that Mi « 0.88, M 2 « 0.42 and M 3 w 0.70 we can conclude that it is 
a genuinely tripartite entangled state. 



4. Maximization of Af-quantities for a three qubit system 

In the second example of section |2] we proved that for the GHZ-state all the M^'s reach 
their maximum value. In this section we prove that all pure states of three qubits for 
which Mfc = 1 Vfc are equivalent to the GHZ-state, i.e. there exists a local and unitary 
transformation that maps this state into GHZ-state. 

Our simple and basic approach makes more transparent the proof of this statement 
in comparison with the proof given in [25] . 
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Consider a three qubits pure state for which it is valid the hypothesis = 1 VA;. 
Since for this state the reach their maximum, we can write the state in the following 
form: 

W = ^|0>*l^>S + ^|l) fc |^) ]t> (17) 

where k refers to the two indexes different from k, {|0) , |1)} is the standard basis of the 
k th qubit and |0 O ) an d \<pi) are such that ((j) a \ (j)b) = 5 a b for a, b = 0, 1, simply applying 
a Schmidt's decomposition. 

By the way, we note that in the case of two qubits, maximization of the unique 
M-functional leads to a Bell-like structure. 

In the tripartite case, because of the simultaneous maximization of the three M^s, 
we have that this structure is valid for each of the three qubits. 

Moreover, this structure is left unchanged by any unitary and local transformation. 

Let us consider this structure with respect to the first qubit and, applying 
again a Schmidt's decomposition, we can convert the state \4>o) 23 m ( cc, s#i |00) 23 + 
e %xi sin^i |11) 23 ), getting: 

M = ±= |0> x (cos*! |00) 23 + e* 1 sm9 1 |11) 23 ) 

+ -Lj* [|0) 2 (a |0} 3 + b |1) 3 ) + |1} 2 (c |0> 8 + d |1) 3 )], (18) 

where {|0);,|1) ; } is a basis, not necessarily standard, for the I th qubit, and a,b,c,d 
are complex numbers and every coefficient satisfies the normalization and orthogonality 
condition between the two states |0o)23 an d |0i) 23 . 

Now we can decompose the state \ip) with respect to the second qubit, getting: 

M = |0} 2 [cos^ |00} 13 + (I), (a |0) s + b |1> 8 )] 

+ -L |1) 2 [e^ sin^ |01) 13 + e*x |l) x (c |0) 3 + d |1> 3 )]. (19) 

The orthogonality conditions involves the orthogonality between the two states 
(a |0) 3 + 6|1) 3 ) and (c|0) 3 + d|l) 3 ), which implies a = acos9 2 , b = ae iX2 sin* 2 , 
c = f3e~ %X2 sin 9 2 and d = (3 cos # 2 . 

In addition, exploiting the normalization constrain with respect to the first qubit, 
we find a = cos6> 3 and (3 = e %xz sin6> 3 , so that we can rewrite the state in equation ([TBI) 
as: 

|^) = i=|O) 1 |0o) 2 3 + -^e iX |l)il0i) 2 3 ) (20) 

where 
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|0o) 23 = (cosfli |00) 23 + e % sintfi |11} 23 ) 

and 



|0i) 23 = cos# 3 |0) 2 (cos# 2 |0) 3 + e lX2 sin# 2 |l) 3 

+ e iX3 sin 6 3 1 1) 2 (e~ iX2 sm6 2 |0) 3 - cosfl 2 |1) 3 ). 

Summarizing, we impose the normalization condition of \<f>o) e with respect to 
the first qubit and the orthogonality condition with respect to the second qubit. We now 
have to enforce: l)The orthogonality condition with respect to the first qubit; 2)The 
normalization condition with respect to the second qubit; 3)The orthogonality condition 
with respect to the third qubit; 4) The normalization condition with respect to the third 
qubit. 

The orthogonality condition for the first qubit leads to e %x cos Q\ cos 9 2 cos # 3 — 
e *(x+X3-xi) gin^! cos# 2 sin# 3 = 0, and so we have the condition: 

02 = | V cos #i cos fl 3 -e i(X3 - Xl) sin^sin^ 3 = 0. (21) 
Rewriting the state \ip) as, 



= -p |0} 2 [costfi |00) 13 + e ix cos# 3 |l> a (cos# 2 |0) 3 + e iX2 sm9 2 |1) 3 )] 
v2 

+ i=|l) 2 [sin^|01) 13 (22) 

+ e i{x+X3) sin 9 3 (e~ lX2 sin9 2 |0) 3 - cosfl 2 |1) 3 )], 

the normalization condition for the second qubit leads to (cos 2 61 + cos 2 9 3 cos 2 9 2 + 
cos 2 6* 3 sin 2 9 2 = 1), that is, 

cos 9 3 = sin 81, (23) 

and likewise for the other state (we consider only this case because possible sign 
difference between sine and cosine can be described by adjusting the phase factors). 
Now, we decompose the state with respect to the third qubit: 

|V> = -^= |0) s l^o>i2 + 4| |l>al#>i2. (24) 

where 



|0' o ) 12 = (cos^i |00) 12 + e ix cos# 3 cos# 2 |10) 
+ e i(x+X3-x 2 ) sin ^ 3 S in 6) 2 |11> 12 ) 
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|0;) l2 = ( e ixi sm6 1 |01) 12 + e i(x+X2) cos^ 3 sin^ 2 |10) 12 
_ e i(x+X3) sin ^ 3 C os 6> 2 1 11> 12 ), 

and we impose the third condition (orthogonality with respect to the third qubit): 
(e* X2 COS3 cos 62 sin 9 2 — e VX2 sin 2 9 3 cos 62 sin #2 = 0), that leads to 

h = \ V 9 2 = V 9 2 = ^, (25) 

(again we consider only 6> 3 = | because possible sign difference between sin 9 3 and cos 9 3 
can be englobed in the phase factors). 

Summarizing, we can have three possible situations: 

V3 — 4 => Vl — 4, 

#2 = 0; (26) 



The normalization condition at point 4, cos 2 9\ + cos 2 9 3 cos 2 #2 + sin 2 #3 sin 2 #2 = 1 
is automatically satisfied for 9 3 = | and #1 = ^ as for # 2 = 0. Instead, if 62 — § then 
necessarily #i = #3 = ~. 

We can apply the conditions we have found, writing the initial state as: 

-L|0) 1 ( s ine 3 |00) 2 3 + e^ cos # 3 |11) 23 ) 

+ -^e ix [cos # 3 |0> 2 (cos # 2 |0) 3 + e iX2 sin# 2 |1) 3 ) (27) 
a/2 

+ e iX3 sin9 3 (e- iX2 sm9 2 |0) 3 - cos# 2 |1> 3 )], 
where parameters vary in the following way: 
03 = f V0 2 ; 

(28) 

#2 = V# 3 ; 

It remains to prove that this state is equivalent to the GHZ under local unitary 
transformation. 

If #2 = 0, condition (J2T1) assures either \i — X3 or $3 = or $3 = f • L e l us examine 
these three possibilities. 

If #2 = and #3 = the state becomes: 

-^=e lX3 1011) + ^=e ix 1100), (29) 
and we can easily find a transformation that sends it in the GHZ. 
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If 2 = and 6*3 = | then we have the state: 
-L|000)-i=e^)| m)) 

and even in this case we can easily find the transformation needed. 
If, finally, 9 2 = and xi — X3 the state becomes: 



-L|0) 1 ( s ine 3 |00) 23 + e ^ cos0 3 |ll> 23 ) 



V2 



e ix \l) 1 (cos 9 3 1 00) 23 - e ixs sin0 3 |11) 23 ) 



and it is converted in the GHZ by: 

sin# 3 |0) 1 + e* x cos# 3 
e ixa cos # 3 10)! - e^+x 3 ) sin# 3 



We have now to examine the first case of eg. (1281) . which is #3 
we have to distinguish the two cases 9 2 = f and 9 2 7^ §. 
If 6 3 — I and 2 7^ f > the state has the form: 



^i»>.(^i«»- 



71' 



11) 



2:! 



+ [_L | ) 2 (cos^ 2 |0) 3 + e^ sinfl 2 |1} 3 ) 

+ -L e *» |1) 2 (e"* 2 sin£ 2 |0) 3 - cos£ 2 |1) 3 )], 



(remember the (}2~Tj) assures xi = X3)- 
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0>i -> 



1) 



1 

V2 



v/2 



0) 2 -> cosf 



0) + e* (xi ~* 2) sin^ 



2 ->■ sin ^ 



0) -e 

2 



0) 3 -> cosf 



1) 3 -> sinf 



1 ) + e 4 * 2 sin f 



0^ -e iXa cosf 



r ) 3 



maps the GHZ-state in the one found. 

If 6*3 = I and 6 2 — f then the state reduces to: 



-L|o)(-L|oo )W *-L|n) 

+ -^e ix II) ( e lX2 ^= 101} + e i(x3 "* 2) 4= 110) 



and we can apply the transformation: 



0>! 



1). 



1 

v/2 



\ + ^ e 



v/2 



1 



v/2 



1 



°} 2 - V' 



J=pi(^-X2 + ^) 



!)a ^2 



0^ + A_S^ + x,-m 

o) 2 - ^^-^ 



that returns the GHZ, and this ends the proof. 
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5. Discussion 



In this paper, by introducing a simple algebraic approach we have reproduced some 
known results which provide recipes to establish whether a certain multipartite pure 
state is entangled with respect to a given bipartition. Such an analysis allows to 
determine whether the pure state describing a system is completely separable, separable 
or totally entangled. In the last case one can say that such a state possesses genuine 
multipartite entanglement. 

Our treatment naturally led us to introduce a class of functionals which include 
quantities traceable back to the concepts of purity and linear entropy. Moreover, 
we have dealt with the problem of the simultaneous maximization of the relevant 
functionals (purities of the reduced matrices), providing an alternative proof to the 
known result that in the case of three qubits the only states that maximize all the 
relevant functionals are the GHZ-state and all the equivalent states (i.e., equal up to 
local unitary transformations). 

This preliminary studies paves the way to the analysis of multipartite entanglement 
in cases wherein more than three subsystems are involved. In particular, simultaneous 
maximization of the relevant functionals for iV-partite systems — with N > 4 — could 
provide interesting results and ideas in the study of maximal multipartite entanglement. 
It is of relevance to stress in addition that other known functionals might be found in 
the class we have introduced, or, alternatively, new quantities could be introduced and 
their properties explored. 
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